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Abstract 

We analyse oscillations observed in the X-ray light 
curve of the late-type star AT Mic. The oscillations oc- 
curred during flare maximum. We interpret these oscilla- 
tions as density perturbations in the flare loop. Applying 
various models derived for the Sun, the loop length and 
the magnetic field of the flare can be estimated. We find 
a period of 740 s, and that the models give similar results 
(within a factor of 2) for the loop length (~ 5.4 10 10 cm) 
and the magnetic field (~ 100 G). For the first time, an 
oscillation of a stellar X-ray flare has been observed and 
results thus obtained for otherwise unobservable physical 
parameters. 
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1. Introduction 

Oscillations in the solar corona have been observed for 
many years. Wavelength regimes ranging from hard X-ray 
right down to radio have been investigated to search for 
evidence of waves. Periods have been found ranging from 
0.02 to 1000 s. Table 1 in |Aschwanden et al. (1999) 



pro- 
vides an excellent summary of the different periods that 
have been found, and an explanation for their existence. 
Most of these waves have been explained by MHD oscil- 



lations in coronal loops. Roberts (2000) has provided an 
excellent review of waves and oscillations in the corona. 

Many of the observations of waves have been deter- 
mined from variations in intensity brightness. Recently, 
however, a huge step forward has been achieved in solar 
coronal physics due to the high spatial resolution avail- 
able with the Transition Region and Coronal Explorer 
(TRACE). The first spatial displacement oscillations have 
been observed in coronal loops l|Aschwan den et al. 1999|l . 
It was suggested that these oscillations were triggered by 
a disturbance from the core flare site. Various MHD waves 
were investigated and it was found that a fast kink mode 
wave provides the best agreement with the observed pe- 
riod of 280±30s. 

One of the most exciting aspects of observing waves 
in this fashion is that it potentially provides us with the 



capability of determining the magnetic field in the corona. 
It is notoriously difficult to measure the magnetic field in 
the corona. Techniques using the near-infrared emission 
lines have been successful, but have poor spatial resolu- 
tion. Most frequently, indirect methods are used such as 
the extrapolations of the coronal magnetic field from the 
photospheric magnetic field that can be measured using 
Zeeman splitting. Nakariakov & Ofman (2001) made use 
of the flare-related spatial oscillations sometimes observed 
when a flare occurs, to determine the magnetic field. They 
assume that the oscillation is due to a global standing kink 
wave, and hence the magnetic field is related to the period 
of the oscillation, the density of the loop, and the length of 
the loop. In the case of |Nakariakov fc Ofman (2001)| they 
found the magnetic field to be 13±9G. The errors on this 
are large due to errors on the determination of density, 
loop length and period. Schrijver & Brown (2000) suggest 
that the oscillations are due to a sudden displacement of 
the magnetic field at the surface which causes an oscilla- 
tory relaxation of the field. Recent work by Nakariakov et al. (2004 
suggests that the oscillations are second standing harmon- 
ics of acoustic waves tied to the loop length. The model 
gives a relationship between the oscillation period, the 
loop length and the average plasma temperature, all of 
which can be independently observed on the Sun. 



Waves are observed across the electromagnetic spec- 
trum on the Sun. However, observations on other stars are 
rare. One reason for this is that on the Sun, we have the 
spatial resolution to zoom in on small regions. For exam- 
ple, the transverse TRACE loops that have been observed 
are away from the main flare site, and hence have lower 
emission measure. Analysing an X-ray light-curve, for ex- 
ample, would be unlikely to provide evidence of waves if 
the emission was averaged over the whole Sun. Oscilla- 
tions have been observed for optical stellar flares (e.g., 



Mathi oudakis et al. 2003l and MullanliTal. 1992jl.|Mathioudakis e 



found a period of 220 s in the decay phase of a white-light 
flare on the RS CVn binary II Peg. Using this they de- 



termine a magnetic field of 1200 G. Mullan et al. (1992) 



again found oscillations in X-ray active red dwarfs. They 
compared their observations to the possibility that they 
were observing p-modes. On the Sun p-modes give rise to 
5-minute oscillations on the surface. They concluded that 
p-modes could not produce the amplitudes they observed 
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and that it must be coronal oscillations, which have, until 
now, not been observed. 

In this work we observe for the first time an oscilla- 
tion in the corona of AT Mic. We measure the periodicity, 
and hence determine the magnetic field. This is the first 
observation of an X-ray oscillation during a stellar flare. 

2. Target 

AT Mic (GJ 799A/B) is an M-type binary dwarf, with 
both stars of the same spectral type (dM4.5e+dM4.5e). 
Both components of the binary flare frequently. The radius 
of AT Mic given by |Lim et al. (1987)|is 2.6 10 10 cm, us - 



ing a stellar distance of 8.14 pc IjGliese fc Jahreiss 1991|) . 
Correcting for the newer value for the distance from HIP- 
PARCOS (10.2±0.5 pc, |Perryman et al. 19971 ), using that 
the stellar radius is proportional to the stellar distance for 
a given luminosity and spectral class, we obtain a stellar 
radius of r* = 3.3 10 10 cm = 0.47 r©. The mass of AT 



Mic given by Lim et al. (1987) is m* = 0.4 toq 



3. Observation 

For our analysis we used the XMM-Newton observations 
of AT Mic on 16 October 2000 during revolution 156. 



Raassen et al. (2003) have analysed this data spectroscop- 
ically, obtaining elemental abundances, temperatures, den- 
sities and emission measures, while a comparative flare 
analysis between X-ray and simultaneously observed ul- 
traviolet emissions can be found in Mitra-Kraev et al. (2004) 
Here, we solely used the 0.2-12 keV X-ray data from the 
pn-European Photon Imaging Camera (EPIC-pn). 

Figure n shows the AT Mic light curve. The observa- 
tion started at 00:42:00 and lasted for 25.1 ks (~ 7 h). 
There is one large flare, starting ~ 4 h into the obser- 
vation, increasing the count-rate from flare onset to flare 
peak by a factor of 1.7, and lasting for 1 h 25 min. It 
shows a steep rise (rise time r r = 1100 s) and decay (de- 
cay time Td = 1700 s). There is an extended peak to this 
flare, which shows clear oscillatory behaviour. The am- 
plitude of the oscillation is around 16%. Applying multi- 



For example, Svestka (1989) shows solar flare light curves 
and a flat top is not seen. An alternative to oscillations for 
this flare is that repeated and rapid flaring is occurring. 
We do not consider this option in this paper, and assume 
for our analysis that we are observing a coronal oscilla- 
tion. The purpose of our analysis is to determine through 
wavelet analysis the period and amplitude of the oscilla- 
tion. We determine the magnetic field and loop length of 
the coronal loop assuming that the oscillation is due to an 
acoustic wave. As a validity check the value of loop length 
was compared to the value determined from a radiative 
cooling model. 



4. Results 
4.1. Oscillation during flare maximum 

To investigate periodicities during the observation, we ap- 
plied a continuous wavelet transformation to the lOs-bin- 
ned-data, using a Morlet wavelet (see, e.g., |Torrence fc Compo 199 
The wavelet coefficients with the lowest periods are dis- 
played in Fig. n] The flare oscillation is picked up with a 
5<7 significance level. The local maximum with the shortest 
period and a significance of more than 5a is simultaneous 
with the flare top and identifies the oscillation, which, at 
the local maximum, has a period of P = 740 s . The period 
interval where the wavelet coefficients have a significance 
of 5a or more, is [730,920] s. 

4.2. Loop length from oscillation 

Interpreting these oscillations as the second spatial har- 
monics of an acoustic wave within the flare IjNakariakov et al. 2004 
we can then determine the flare loop length from 



L 



1 



[Mm] « — • P [s] • VTp] . 



(1) 



temperature fitting, Raassen et al. (2003) obtain a mean 
flare temperature T = 24 ± 4 MK and a quiescent tem- 
perature T q = 13 ± 1 MK, and from the O vn line ratio a 
flare and quiescent electron density of n = 4+5 10 10 cm -3 



and n a 



1.9 ± 1.5 10 10 cm 



respectively. The total 



flare and quiescent emission measures are EM = (19.5 ± 
0.8)10 51 cm" 3 and EM q = (12.2 ± 0.5)10 51 cm" 3 . 

One of the reasons why we concentrated on this obser- 
vation is that following the rise phase of the flare there was 
a period of approximately 40 minutes when the light curve 
stayed at a high intensity level and showed strong oscilla- 
tions. The oscillatory activity can be clearly seen by eye 
(see Fig.QJ, and shows behaviour suggestive of damping. 
This character is much different from solar flares where 
the flare reaches a peak rapidly, followed by a slow decay. 



Inserting the above obtained period P — 740 s and a 
flare temperature T — 24 MK, we obtain a loop length 
L = 5.4 10 10 cm. Note that the temperature carries a la 
error of ~ 16%, therefore the error of L is at least 8% 
(AL > 0.5 10 10 cm). It is not straightforward to calculate 
the corresponding la error of P. With these numbers, the 
inferred speed of sound is c s = L/P w 730 km/s. 

4.3. Loop length from radiative cooling times 

The loop length can also (and independently of any os- 
cillation) be estimated from rising and cooling times ob- 
tained from the temporal shape of the flare, applying a 
flare heating/cooling model (see, e.g.,|Cargill et al. 19951 ). 
We follow the approach by Hawley et al. (1995)| who~m- 
vestigated a flare on AD Leo observed in the extreme ul- 
traviolet. The shape of this flare is very similar to our 
flare on AT Mic, but roughly 10 times larger (in duration 
as well as in the increase of the count rate). It also shows 
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Figure 1. The light curve and wavelet coefficients of AT Mic. The observation started on 2000-10-16 at 00:42:00 h and 
lasted for 25'100 s (~ 7 h). The upper panel shows the observed light curve, binned up to 100 s. The vertical dashed 
lines indicate the start and end of the rise and decay phase of the flare. The lower panel shows the absolute values of a 
section of the corresponding wavelet coefficients (see main text) divided by their standard deviation. The 3, 4 and 5 a 
significance contour lines are drawn. The two dashed lines mark the border of the cone of influence. A local maximum 
is clearly seen with a period of around 740 s during the flare peak. 



a flat top with a possible oscillation. The flare loop energy 
equation for the spatial average is given by 

\i> = Q-R, (2) 

with Q the volumetric flare heating rate, R the optically 
thin cooling rate and p the time rate change of the loop 
pressure. During the rise phase, strong evaporative heating 
is dominant (Q > JJ), while the decay phase is dominated 
by radiative cooling and strong condensation (R ^> Q). At 
the loop top, there is an equilibrium (R = Q). The loop 
length can be derived as 



L = 



1500 



4/7 



„1.58A 4 / 7 d 



r 3/7 . T l/2 



(3) 



where r r is the rise time, Td the flare decay time (indicated 
in Fig. with the vertical dotted lines), T the apex flare 
temperature and x 2 d = Cd/c maxi with c max the peak count 



rate and Cd the count rate at the end of the flare. Inserting 
these values, the loop length becomes L = 2.3 10 10 cm, 
which is about a factor of 2 smaller that the loop length 
inferred from the oscillation. 



4.4. Determination of the Magnetic field 
Interpreting the oscillations in terms of global standing 



kink waves, Nakariakov & Ofman (2001) derive a relation 
for the magnetic field 



B 




(4) 



with p the mass density inside, and p ext the mass density 
outside the loop. Note that there is a correction factor of 



0.64 in the equation as presented by Nakariakov & Ofman (2001) 
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UMathioudakis et al. 2003f) . We can simplify Eq. |@} by in- 
serting Eq. and using p = m p ■ n and p ext — m p ■ n q , 
with m p the proton mass. Then, the magnetic field is given 
by 

B = #- 10S '\/ T "V>( 1 + ^)' (5) 

Inserting the temperature and densities from the spec- 
troscopy results, we obtain a magnetic field of B = 115 ± 
56 G. Note that the large error is resulting mainly from the 
density uncertainty. Using the loop length from the radia- 
tive cooling approach fSect. |4~S|) and Eq. the magnetic 
field is B « 50 G. 

4.5. Pressure balance 

To maintain stable flare loops, the gas pressure of the 
evaporated plasma must be smaller than the magnetic 
pressure 

B 2 

2nkT < — . (6) 

Knowing the flare density and temperature, we get a lower 
limit for the magnetic field B > 80±60 G. Again, there is a 
large error because of the large uncertainty for the density. 
Shibata &: Yokoyama (2002)| assume pressure balance and 
give equations for B(EM,n q ,T) and L(EM,n q ,T) (their 
Eqs 7a, b). Using these relations, we obtain a magnetic 
field of B = 70 ± 40 G and a loop length of L = (2.8 ± 
1.7) 10 10 cm. 

5. Discussion and Conclusions 

We have used three different approaches to determine the 
loop length of the AT Mic flare and two different ap- 
proaches for the magnetic field. The loop length derived 
in Sect. 14.21 from the flare oscillation is the largest with 
L = 5.4 • 10 10 cm, while the loop length derived from ra- 
diative cooling times (Sect. I4.3fl is somewhat smaller by 
roughly a factor of 2, L = 2.3 10 10 cm. The loop length 
derived from pressure balance fSect. l4~5*| is similar to the 
latter one, L = (2.8±1.7) 10 10 cm. The loop length derived 
by assuming that the variations seen are due to oscillations 
is a factor of 2 different to the loop length determined by 
more usual methods. Considering the assumptions that 
have to be made, this is good agreement, and gives us 
confidence that we are, for the first time, observing a stel- 
lar coronal loop oscillating. The magnetic field ranges from 
50 G (radiative cooling) to 115±56 G (oscillation), with 
the value from pressure balance in between with 70±40 G. 
Again, these three values are all consistent. 

This was the first time that an oscillation during flare 
peak was observed in X-rays in a stellar flare and flare 
loop length and magnetic field derived from it. The values 
are consistent with other flare models. 

Acknowledgements 
We acknowledge financial support from the UK Particle Phy- 
sics and Astronomy Research Council (PPARC). UMK would 



also like to thank the European Space Agency (ESA) and the 
University College London (UCL) Graduate School for finan- 
cial assistance to attend the Cool Stars 13 conference. 

References 

Aschwanden M.J., Fletcher L., C. J. Schrijver et al. 1999, ApJ 
520, 880 

Cargill P.J., Mariska J.T., Antiochos S.K. 1995, ApJ 439, 1034 
Gliese W., Jahreiss H. 1991, Preliminary Version of the Third 
Catalogue of Nearby Stars, Astron. Rechen-Institut, Hei- 
delberg 

Hawley S.L., Fisher G.H., Simon T. et al. 1995, ApJ 453, 464 
Lim J., Nelson G.J., Vaughan A.E. 1987, Proc. ASA 7, 2 
Mathioudakis M., Seiradakis J.H., Williams D.R. et al. 2003, 
A&A 403, 1101 

Mitra-Kraev U, Harra L.K., Gudel M. et al. 2004, A&A, sub- 
mitted 

Mullan D.J., Herr R.B., Bhattacharyya S. 1992, ApJ 391, 265 
Nakariakov V.M., Ofman L. 2001, A&A 372, L53 
Nakariakov V.M., Tsiklauri D., Kelly A. et al. 2004, A&A 414, 
L25 

Perryman M.A.C., Lindegren L., Kovalesky J. et al. 1997, A&A 
323, L49 

Raassen A.J.J. , Mewe R., Audard M. et al. 2003, A&A 411, 
509 

Roberts B. 2000, Sol. Phys. 193, 139 
Schrijver C.J., Brown D.S. 2000, ApJ 537, L69 
Shibata K., Yokoyama T. 2002, ApJ 577, 422 
Svestka, Z., 1989, Sol. Phys. 121, 399. 

Torrence C, Compo G.P. 1998, Bull. Amer. Meteor. Soc. 79, 
61 



